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Приведена постановка и схема реализации нового метода решения плоских задач теории упругости, основанного на диагона- 
лизации системы уравнений равновесия. Получены аналитические решения трех плоских задач в напряжениях для нагружения 
полосы сложной нагрузкой. 


Введение 

При решении двумерных задач теории упруго- 
сти 11-4], возникает проблема отыскания двух не- 
известных функций (перемещений и, ѵ в направле- 
нии координатных осей ОХ, ОГкоторые удовлетво- 
ряют системе двух дифференциальных уравнений в 
частных производных второго порядка - системе 
уравнений равновесия в форме Ляме [1]. Кроме то- 
го, искомые перемещения должны удовлетворять 
заданным условиям на граничном контуре. В клас- 
сической постановке задачи теории упругости в на- 
пряжениях граничные условия относительно иско- 
мых перемещений так же, как и уравнения равно- 
весия, образуют систему двух дифференциальных 
уравнений в частных производных (в отличие от 
системы уравнений равновесия система граничных 
условий имеет первый порядок). 

Для решения плоских задач теории упругости 
предлагается использовать метод, основанный на 
преобразовании системы уравнений равновесия к 
диагональному виду [5, 6]. Очевидно, что при та- 
ком подходе к отысканию решения появляются не- 
которые преимущества, облегчающие процедуру 
интегрирования уравнений равновесия и удовле- 
творения граничным условиям. Данные преиму- 
щества можно использовать как в процессе постро- 
ения аналитических решений, так и при численной 
реализации с помощью конечно-разностного, ко- 
нечно-элементного и других методов. В настоящей 
работе применяется аналитический метод; резуль- 
таты сопоставляются с известными решениями, 
полученными с помощью функции напряжений. 

1. Постановка плоской задачи теории упругости 
для диагонализованной системы уравнений 
равновесия и схема решения 


АѲ = 0, (1) 

Ак = 0. (2) 

Здесь Ѳ - объемная деформация, которая выра- 
жается через перемещения согласно 

Ѳ = б х и + <і 2 ѵ. 
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где к - произвольная гармоническая функция, А - 
оператор Лапласа 

А = б 1 2 + 6 2 , 
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б) Граничные условия. 

Для пары новых искомых функций Ѳ=Ѳ(х,у), 
к=к(х,у) граничные условия преобразуются в усло- 
вия вида: 

а х 1 + т ху т = х п , 

а у т + т^ = У п , (3) 

где Х„, Ѵ„ - компоненты объемных сил, сг г , а у , т ч - 
напряжения; I, т - направляющие косинусы нор- 
мали к граничному контуру, имеющему дугу .у: 

/ = соз(п,х) = — , т = софт, у) = - — . 
сіз 

в) Функции, гармонически-сопряженные . 

Кроме Ѳ=Ѳ(х,у) и к=к(х,у), для расчета компо- 
нент напряженно-деформированного состояния 
необходимо задание еще одной пары функций, 
удовлетворяющих системе Коши-Римана 

сі х Ѳ = асі 2 со, (1 2 Ѳ = -асіуо. ( 4 ) 


Рассмотрим основные элементы предлагаемой 
модификации постановки плоской задачи, осно- 
ванной на диагонализации системы уравнений 
равновесия в форме Ляме. 

а) Уравнения равновесия в области 
Вместо отыскания пары функций и=и(х,у), 
ѵ=ѵ(х,у), удовлетворяющих системе уравнений рав- 
новесия Ляме, модифицированная постановка 
предполагает решение двух гармонических уравне- 
ний: 


Здесь а = — — , со — компонента вектора вращения 
1 + А 


Исходя из соотношений (4), функцию асо{х,у), 
гармонически-сопряженную Ѳ(х,у), можно опреде- 
лить в любой точке расчетной области согласно 

(Х,У) 

асо(х,у)= | (-с/ 2 Ѳ(х,у)с/х + с/ І Ѳ(х,у)с/у) + С„ (5) 
0 
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где С ( - константа. 

Аналогично можно найти функцию %=%(х,у), 
гармонически-сопряженную к(х,у): 

<і 1 к = -<і 2 х, й г к=й х х, (6) 

и, таким образом, х( х >У) в любой точке расчетной 
области: 

и,у) 

Х(х,у) = | (<і 2 к(іх-(і 1 к<іу)+С 2 , (7) 

0 

где С 2 - константа. Далее будем использовать зна- 
чения С|,С, равные нулю. 

г) Формулы для компонент напряженно-деформи- 
рованного состояния. 

По найденным из решения двух краевых задач 
типа Дирихле функциям Ѳ, к и по определению 
функций асо и х, гармонически-сопряженных Ѳ, к, 
можно рассчитать напряжения в расчетной области: 

<7 Г =Ѳ-^[хД,Ѳ -ус1 2 Ѳ\ + к, 

<т ѵ , =Ѳ + ^[хс/,Ѳ -ус1 2 Ѳ]-к, 

х ху =-\[хсі 2 Ѳ + ус1 і Ѳ] + х- ( 8 ) 

Перемещения и, ѵ упругого тела можно рассчи- 
тать по напряжениям (8), пользуясь стандартными 
соотношениями [2]. 

Таким образом, реализация модифицированной 
постановки плоской задачи теории упругости, 
определяемая соотношениями (1-8), имеет некото- 


рые отличия от классической постановки. Основ- 
ное преимущество предлагаемой модификации за- 
ключаются в том, что уравнения равновесия имеют 
диагональный вид (1,2). Действительно, уравнения 
(1,2) представляют собой дифференциальные ура- 
внения в частных производных второго порядка, а 
не систему уравнений равновесия в форме Ляме. 

Для иллюстрации метода рассмотрим следую- 
щие задачи. 

2. Нагружение пластины сложной нагрузкой 

Пластина размером И, Ь подвергается действию 
напряжений, эпюры которых представлены на рис. 1. 

Для решения необходимо записать аналитиче- 
ские выражения граничных напряжений. Сначала 
получим вид функции У„(х,у). 

а) Рассмотрим граничные значения а; на верхней 
и нижней гранях, обозначив их через а 

Н а 1 ,3 „ (ИЛ* 

при у = — а = —ап =2 а\ — \ , т = 1; 

2 у 4 \2 ) 

И а 1 ,3 „ (НЛ 3 

при у = — <т „= — ап = -2а — , т = -1. 

2 ѵ 4 [2] 

Здесь й=соп8і. 

Объединяя два последних выражения, получим 

сг“=2 ау 3 т\ н . 

у І,=±- 

б) Рассмотрим вклад У„(х,у) в касательных напря- 
жений на правом торце. Согласно эпюре гра- 
ничного т„ имеем: 




Рис 1. Эпюры нагружения пластины сложной нагрузкой 
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при у = — г, = — аЫг =-6аЬ 
2 " 2 

А 3 2 

при у = — г = — аЬп =-ЬаЬ 
2 ™ 2 

при у = 0 = О. 

Значит, зависимость т ху от по торцу выражается: 
т ху =-6аЬу 2 . 

Далее учтем изменение от левого торца к пра- 
вому. Поскольку на левом торце (х=0, /=-1) г ч =0, а 
на правом (х=Ъ, 1=1), т ч =—6 аЪу 1 , то получаем: 

т ху =-6аху 2 . 

Таким образом, функция Ѵ„ имеет вид: 

Т =а(2у г т-6ху 2 1). 

Аналогичным образом рассчитывается анали- 
тическое выражение для Х„: 

Х„(х,у) = а[(6х 2 у-4у 3 )1-6ху 2 т]. 
Зададимся функциями Ѳ(х,у ) и к(х,у) вида: 

Ѳ(х,у) = А(Зх 2 у-у 3 ), (9) 

к(х,у) = В(Зх 2 у-у 3 ), (10) 

где А и В - неизвестные константы, которые 

необходимо найти. Функции (9, 10) соответствен- 
но удовлетворяют условиям (1, 2), то есть являются 
гармоническими . 

Функция %(х,у) определяется из (10) с помощью 
(7) (для упрощения записи здесь и далее, коорди- 
наты х,у будем опускать): 

Х = В(х 3 -Зху 2 ). (11) 

Данная функция тоже является гармонической. 
С учетом (9-11) распишем напряжения (8): 

а у =[^А-ЗВуу + [^А + В^у 3 , 

1 х? =^--А-ЗВ^ х у 2 +^--А + В^х г . ( 12 ) 

Рассмотрим левую границу пластины, х=0. 
Здесь /= — 1, т= 0, соответственно 

х п = 4 а у 3 , 

У„= 0 , ( 13 ) 

а граничные условия (3) преобразуются так: 

~Ѵ х = Х п’ 

^ ху = У. (14) 

Подставляя в (14) функции из (12) и (13) и ре- 
шая полученные уравнения, получаем, что 

В = 4а-^А. (15) 




Учитывая (15), напряжения (12) будут иметь 
следующий вид: 

а х = (-6А + 12а)х 2 у-4а у 3 , 
а у = (12А-12а)х 2 у+ (-2А+ 4а)у 3 , 
т ху = (-9 А - 12 а)ху 2 + (-4А+ 4 а)х 3 . (16) 

Теперь рассмотрим правую границу пластины, 
х=Ъ. Здесь 1=1, т=0, соответственно 

Х п =а(6Ь 2 -4у 3 ), 

У п =-6аЬу 2 , (17) 

а граничные условия (3) будут выглядеть так: 

=*,. 

Х у = У п - ( 18 ) 

Подставляя в (18) функции из (16) и (17) и ре- 
шая полученные уравнения, получаем, что 

А = а, (19) 

далее, по ( 15) с учётом (19) определяем, что 

В = \а. ( 20 ) 

Таким образом, обе неизвестные константы 
найдены. Расписывая (9, 10) с учётом (19, 20), по- 
лучаем, что Ѳик имеют вид: 

Ѳ = а{Зх 2 у-у 3 ), 

к = ^а(Зх 2 у-у 3 ). (21) 

Теперь определены все функции, необходимые 
для нахождения напряжений. Подставляя функции 
(21) в (8), получаем: 

сг х (х,у)=6ах 2 у-4ау ъ , 
а у (х,у) = 2 а у 3 , 

т ху (х,у) = -6аху 2 . (22) 


Сравним полученное решение с классическим, 
которое можно рассчитать с помощью функции 
напряжений. Для данной задачи функция напря- 
жений выглядит следующим образом [4]: 


ср = а 


х 2 у 3 - 


_У_ 

5 


соответственно, напряжения будут следующие: 
а х = с1 2 2 (р = 6ах 2 у-4ау 3 , 

(У у = й 2 ц> =2ау ъ , 

т ху =-<і і (і 1 (р = -ваху 2 . (23) 


Видно, что соответствующие функции в (22) и в 
(23) попарно совпадают. 


3. Нагружение пластины, силами, 
распределенными по линейному закону 

Пластина размером И, Ъ подвергается действию 
напряжений, эпюры которых представлены на 
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рис. 2. Аналогично предыдущей задаче определяет- 
ся, что нагрузка в функции координат имеет вид: 
х „ О, У) = 0, У„ (х,у) = вахт. 


Объединяя (27) с (29) и (28) с (30), получаем, что 
о = -\в,с = - х -а, (31) 


ік У 

^^тттГТТТТПШШ 

+6 аЬ 

А А 4 1 1 

А V 

0 

А V 
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уг 

'**І1 
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Рис. 2. Эпюры нагружения пластины силами, распределен- 
ными по линейному закону 


Для решения данной задачи зададимся гармо- 
ническими функциями Ѳи (св виде: 

Ѳ = Ах + Ву, 

к=Сх + Оу, (24) 

где А, В, С, В - константы, которые будет необхо- 
димо найти. По (7) находим, что 

X = В)х-Су. (25) 

Все три функции, как легко показать, являются 
гармоническими . 

Используя (8) получаем следующий вид для на- 
пряжений: 

(т х = х[^А + С^ + у^В + Б^, 

т *>= х [~\ В + 0 ] + у [-\ А ~ С ]- (26) 

Рассмотрим левую границу пластины. Для этой 
границы х=0, /=- 1 , т={). С учётом этого, из (3) по- 
лучаем: 

°х = 0, (27) 

^. = 0, (28) 

при этом, учитывая (26) на рассматриваемой гра- 
нице, с7 д . и г х>! принимают вид: 

а х =у\^В + 0^, (29) 

Т ^ = У І~2 А ~ С \ (30) 


и, с учётом этого, функции из (26) принимают сле- 
дующий вид 

ст д = 0, 

су у = 2Ах + 2В у, 

Т ху = ~2В х. (32) 

Теперь рассмотрим верхнюю границу пласти- 
ны. Для этой границы У = — , / = 0, т = 1, а функ- 
ции (3) имеют следующий вид: 

о у = бах, (33) 


= 0, 


при этом, учитывая (32), 


а =2Ах + 2В — , 
2 


(34) 

(35) 


т ху = -2Вх. (36) 


Объединяя (34) и (36), получаем, что 5=0, и, 
учитывая это при объединении (33) и (35), получа- 
ем, что А=3а. 


Из (31) теперь находим, что С = -— а,0 = 0. 


Подставляя значения найденных констант в 
(32), находим, что 

С7 Г = 0, 

<т ѵ =6 ах, 

т ѵ =0. (37) 


Сравним полученное решение с классическим. 
В задаче для треугольной нагрузки бигармониче- 
ская функция напряжений имеет вид ср=а> ? соот- 
ветственно, напряжения будут следующие: 

а х (х,у) = с( 2 ц> = В, 2 (ах 2 ) = 0, 

о Г (х, у) = В 2 ср = с1 2 (ах ъ ) = <7,(3ах 2 ) = 6 ах, 

т (х,у) = -В х В,ср = -В х В,(ах 2 ) = 0. (38) 


Таким образом, напряжения, рассчитанные по 
(37, 38), совпадают между собой. 


4. Расчёт плотины треугольного профиля 

Плотина, подвергающаяся нагрузке от давле- 
ния воды, имеет треугольное сечение, как показано 
на рис. 3. Способом, использованным в предыду- 
щих двух задачах, определяется, что нагрузка в 
функции координат имеет вид: 

Х „(х,у) = у уІ + у\% 2 (а)хт, 

У„ (х, у) = (-2/ 1§ 2 (а)х- у 1§ 2 (а) у) т + у 1§ 2 (а)х I. 

Здесь й=соші, у=сопм. 
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Как и в предыдущей задаче, зададимся гармо- 
ническими функциями Ѳм къ виде (24). Функция 
X в этом случае будет иметь вид (25), а напряжения 
- (26). Как и в предыдущем случае, необходимо 
найти значения констант А, В, С, О. 



Рис. 3 . Эпюры нагружения плотины 


Рассмотрим вертикальную поверхность (левую 
границу) плотины. Для этой границы х=0, /=— 1, 
т= 0. Из (3) определяем, что 

о-, = ГУ, (39) 

^=0, (40) 

при этом, <7 Х и т ѵ на рассматриваемой границе, 
принимают вид (29, 30). Объединяя попарно (39) с 
(29) и (40) с (30), получаем, что 

0 = у~В,С = ~А, (41) 

и, с учётом этого, функции из (26) принимают сле- 
дующий вид 

= УУ> 

(Гу =2Ах + у(2В-у), 

т ху =х{-2 В + у). (42) 

Теперь рассмотрим наклонную поверхность 
(правую границу) плотины. Для этой границы рас- 
пишем: х=-ус{%(а), /=ып(а), /я=со8(а). Компо- 
ненты объемных сил Х„ и У„ здесь равны нулю. С 
учетом этого, граничные условия в форме Коши 
будут иметь следующий вид: 
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ст х 8Іп(а) + т соз(а) = О, 

<т со8(а) + г ѵѵ 8Іп(а) = 0. (43) 

Рассмотрим первое уравнение системы (43). 
Подставляя в него функции из (42), получаем: 
у^8Іп(а) -(-25 + у)усІ§(а)со8(а) = О, 

откуда получаем, что 

В = - Г -(^\а)~ 1). (44) 

Далее рассмотрим второе уравнение системы 
(43). Подставляя в него функции из (42), получаем: 
-2 Ау с(§(а ) со8(а ) + у(2В - у) сов (а) - 
-у(-2В + у ) сіё(а ) 8іп(а) = О, 
откуда с учетом (44) получаем, что 

А = -уів 3 (а). (45) 

По (41) и (44) определяем, что 
В=- У -(іё\а)~ 1), 

с=и%\а\ 

0 = У ~( 1 + 31§ 2 (а)). (46) 

Подставляя (45) и (46) в (26), получаем: 

=УУ> 

<у у =-2хущ\а)-уу1ё 2 (а), 

^ Х у=хуіё 2 (ос). 

Сравнение данных результатов с классически- 
ми показало их полное совпадение. 

Вывод 

Показано, что в отличие от классической поста- 
новки плоской задачи теории упругости в случае 
заданных на границе напряжений, модифициро- 
ванная постановка позволяет разделить (для каж- 
дой из двух искомых функций) процедуру опреде- 
ления решения, удовлетворяющего уравнению 
равновесия в области и граничным условиям типа 
Коши на контуре. Метод проиллюстрирован на 
аналитических решениях плоских задач. Результа- 
ты показали полное совпадение с решениями, по- 
лученными с помощью классического метода. 
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